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I.     INTRODUCTION 

A.         BACKGROUND 

Electronic  equipment  has  made  its  way  into  practically  every  aspect  of  modern  life, 
from  high-powered  computers  to  toys  and  appliances.  The  reliability  of  the  electronics  of 
a  system  is  a  major  factor  in  the  overall  reliability  of  the  system.  Electronic  components 
become  potential  sites  for  excessive  heating  due  to  the  fact  that  they  depend  on  passage  of 
electric  current  to  operate.  Continued  miniaturization  of  electronic  systems  has  resulted  in 
a  dramatic  increase  in  the  amount  of  heat  generated  per  unit  volume.  As  stated  by  Cengel 
[Ref.  1],  these  heat  fluxes  are  comparable  in  magnitude  to  those  encountered  at  nuclear 
reactors  and  the  surface  of  the  sun.  These  high  rates  of  heat  generation  can  result  in  high 
operating  temperatures  (in  excess  of  80°  C)  for  electronic  equipment,  which  could 
jeopardize  its  safety  and  reliability.  This  excessive  heat  generation  must  be  properly 
designed  for  and  controlled  to  ensure  the  reliability  of  the  electronic  equipment. 

Natural  and  forced  convection  cooling  with  air  has  been  the  mainstay  method  for 
providing  the  necessary  cooling  of  electronic  equipment.  The  relatively  lower  heat  transfer 
capabilities  of  these  methods  has  led  designers  to  investigate  the  use  of  immersion  cooling 
with  liquids.  With  immersion  cooling  the  heat  generating  components  are  totally  immersed 
in  a  liquid  which  is  the  cooling  medium.  Immersion  cooling  is  better  able  to  withstand  the 
high  heat  flux  and  volumetric  heat  generation,  which  are  characteristic  of  high  performance 
computers.  Immersion  cooling  provides  two  solid  advantages  over  natural  and  forced 
convection  cooling  with  air:  first,  liquids  as  coolants  are  more  effective  than  gases  due  to 
higher  heat  capacities;  secondly,  cooling  with  liquid  has  the  additional  possible  capability  of 


two  phase  boiling  heat  transfer.  Large  amounts  of  heat  can  be  extracted  due  to  the  latent 
heat  of  boiling. 

If  the  problem  of  cooling  an  electrically  heated  wire  immersed  in  an  otherwise 
quiescent  pool  of  liquid  at  its  saturation  temperature  is  examined,  several  heat  transfer 
mechanisms  are  encountered  as  the  wire  temperature  is  increased.  When  the  wire 
temperature  is  only  slightly  higher  than  the  liquid  saturation  temperature  the  only  heat 
transfer  mechanism  is  natural  convection  in  which  buoyant  forces  cause  the  liquid  to  rise 
and  cool  the  wire.  This  is  a  relatively  inefficient  cooling  method.  This  process  has  been 
extensively  studied  and  a  correlation  has  been  established  to  predict  the  amount  of  heat 
transferred  due  to  natural  convection.  Kuehn  and  Goldstein  [Ref  3]  established  the  most 
recent  natural  convection  heat  transfer  correlation  for  small  diameter  horizontal  circular 
cylinders. 

Dhir  states  [Ref.  2],  as  wall  superheat  (difference  between  heated  surface 
temperature  and  liquid  saturation  temperature;  (TSUIf  -  Tsat))  increases,  larger  vapor  bubbles 
begin  to  appear  on  the  heated  surface.  This  is  referred  to  as  the  onset  of  nucleate  boiling. 
Vapor  or  gas  trapped  in  imperfections  such  as  cavities  and  scratches  on  the  heated  surface 
serve  as  nucleation  sites  for  bubbles.  A  wedge  shaped  imperfection  on  a  surface  will  trap 
vapor  or  gas  as  long  as  the  static  contact  angle  between  the  solid  surface  and  the  liquid- 
vapor  interface  is  greater  than  the  wedge  angle.  A  bubble  will  form  if  the  temperature  of 
the  liquid  at  the  tip  of  the  cavity  is  at  least  equal  to  the  saturation  temperature 
corresponding  to  the  vapor  pressure  in  the  bubble.  The  balance  of  forces  that  act  on  the 
bubble  dictates  the  diameter  to  which  the  bubble  will  grow  before  departing.  These  forces 


are  associated  with  the  inertia  of  the  liquid  and  vapor,  the  liquid  drag  on  the  bubble, 
buoyancy,  and  the  surface  tension.  When  the  buoyancy  of  the  bubble  exceeds  the  surface 
tension  of  the  liquid  and  liquid  drag  forces,  the  bubble  is  able  to  depart  the  heated  surface 
and  rise  in  the  fluid.  After  departure,  cooler  liquid  from  the  bulk  fills  the  space  vacated  by 
the  bubble.  When  the  required  superheat  is  once  again  achieved  at  the  tip  of  the  cavity,  a 
new  bubble  starts  to  form  at  the  same  nucleation  site  and  the  bubble  growth  and  departure 
sequence  repeats.  This  repeating  process  is  nucleate  boiling. 

For  highly  wetting  liquids,  only  sites  of  very  small  size  will  be  active  in  bubble 
generation.  For  these  liquids,  the  expected  wall  superheat  at  nucleation  should  be  higher 
than  for  partially  wetting  liquids.  This  excess  superheat  necessary  to  establish  nucleate 
boiling  creates  a  hysteresis  in  the  boiling  curve  since  the  heated  surface  temperature  drops 
once  nucleate  boiling  is  fully  developed.  This  phenomenon  only  occurs  during  the 
transition  from  natural  convection  heat  transfer  to  nucleate  boiling  heat  transfer.  During 
the  transition  from  nucleate  boiling  to  natural  convection  heat  transfer,  there  is  no 
temperature  overshoot  and  the  heated  surface  temperature  steadily  decreases  with 
decreasing  wall  heat  flux.  The  required  excess  superheat  necessary,  to  commence  nucleate 
boiling  is  the  major  problem  with  using  this  method  for  electronic  cooling  applications. 

Immersion  cooling  research  using  dielectric  liquids  such  as  Fluorinerts  (FC  series  of 
3M  Corporation)  is  still  ongoing.  The  dielectric  liquids  are  a  synthetically  manufactured 
fluid  derived  from  common  organic  compounds  by  replacement  of  all  carbon-bound 
hydrogen  atoms  with  fluorine  atoms.  [Ref.  4]  The  dielectric  liquid  used  in  numerous 
boiling   studies   is  FC-72.      FC-72   has   excellent   electrical   insulating   and   dielectric 


characteristics.  The  dielectric  strength  for  liquid  FC-72  is  38kV/0.1  inch  and  the  vapor 
phase  has  a  similar  value.  FC-72  is  non-corrosive  and  has  a  relatively  low  saturation 
temperature  of  56  °C  at  atmospheric  pressure.  [Ref.  4]  These  properties  of  FC-72  make  it 
an  excellent  choice  for  electronics  cooling. 

Conducting  boiling  experiments  with  FC-72  does  have  some  problems.  FC-72  has 
very  low  surface  tension,  which  means  it  is  highly  wetting,  and  therefore  has  the  ability  to 
flow  into  every  small  crack  or  crevice  in  the  heated  surface.  This  makes  the  initial 
formation  of  bubbles  very  difficult  and  will  result  in  a  large  wall  superheat  to  initiate 
nucleate  boiling. 

The  present  study  is  concerned  with  investigating  methods  to  achieve  a  reduction  in 
the  wall  superheat  required  to  initiate  nucleate  boiling.  The  wall  superheat,  also  called 
temperature  overshoot  (TSUrf  -  Tsat),  is  a  major  drawback  to  using  immersion  cooling  with 
electronic  equipment.  The  electronic  equipment  can  not  withstand  the  high  temperatures 
created  prior  to  nucleate  boiling.  One  method  used  to  enhance  the  onset  of  nucleate 
boiling,  which  will  also  reduce  the  temperature  overshoot,  is  the  mechanical  removal  of  the 
bubbles  from  the  heated  surface.  Removal  of  the  bubbles  prior  to  their  normal  departure 
will  create  more  nuclei  sites  at  a  lower  wall  temperature.  This  action  will  facilitate  the 
onset  of  nucleate  boiling  and  create  a  lower  required  temperature  overshoot. 
B.         PREVIOUS  WORK 

You,  Bar-Cohen,  and  Simon  [Ref.  5]  conducted  experiments  with  horizontal 
cylindrical  heater  surfaces  immersed  in  FC-72  and  R-113,  saturated  at  one  atmosphere 
pressure.  They  observed  the  effects  of  fluid  and  surface  properties  and  materials  on  boiling 


regimes  and  on  the  onset  of  nucleate  boiling.  They  found  the  incipience  wall  superheat  for 
FC-72  was  about  25  °C  lower  and  nucleate  boiling  superheat  about  6  °C  lower  than  those 
values  for  R-l  13.  They  have  concluded  FC-72  is  more  advantageous  for  cooling  than  R- 
113. 

You,  Simon,  and  Bar-Cohen  [Ref  6]  studied  enhanced  surfaces  with  more 
available  nucleation  sites  in  immersion  cooling  with  boiling.  They  were  attempting  to 
modify  the  surface  of  the  computer  chip  to  enhance  the  onset  of  nucleate  boiling.  A 
particular  layering  technique  was  found  to  lower  the  wall  superheat  in  boiling  when 
compared  to  the  untreated  surface.  The  new  particle  layering  technique  could  be  applied 
to  computer  chips  without  inducing  stress  or  damaging  the  chip. 

You,  Ammermam,  and  Hong  [Ref  7]  conducted  a  series  of  pool  boiling 
experiments  on  saturated,  gas  saturated,  and  pure  sub-cooled  FC-72.  The  dissolved  gases 
significantly  enhanced  the  transition  to  nucleate  boiling  in  FC-72.  For  the  present  study, 
You's  findings  demonstrate  the  need  to  ensure  that  FC-72  has  been  entirely  degassed  prior 
to  starting  the  experiment.  This  is  critical  since  the  only  nucleate  boiling  enhancement 
desired  is  fluid  oscillation,  which  will  remove  the  bubbles  from  the  platinum  wire.  Any 
remaining  dissolved  gases  in  the  FC-72  will  affect  the  nucleate  boiling  and  skewer  the 
findings  of  this  study. 

Marcellus  [Ref.  8]  conducted  a  series  of  experiments  with  FC-75  to  prove  that 
forced  flow  can  increase  the  critical  heat  flux  value.  This  can  be  accomplished  by  the 
rapid  removal  of  the  vapor  bubbles  from  the  heated  surface. 


Kelleher,  et  al.  [Ref.  9]  conducted  experiments  on  a  vertical  array  of  five 
horizontal  platinum  wires  mounted  in  a  pool  of  FC-72.  Excess  wall  superheat  was 
necessary  to  initiate  nucleate  boiling  in  the  FC-72.  Their  results  show  that  the  bubbles 
from  the  upstream  source  reduced  the  temperature  overshoot  on  the  downstream  surface. 
Using  the  wake  flow  from  one  heated  surface  to  assist  in  removing  bubbles  from  another 
heated  surface  proved  to  be  a  successful  method  for  lowering  the  amount  of  wall 
superheat  necessary  to  initiate  nucleate  boiling. 

Turk  [Ref.  10]  studied  the  effects  of  oscillating  the  fluid  on  reducing  the  excess 
superheat.  Turk  conducted  a  series  of  experiments  using  a  platinum  wire  immersed  in  a 
main  chamber  filled  with  FC-72.  A  smaller  front  chamber  was  rigidly  mounted  to  the  main 
chamber.  The  front  chamber  and  main  chamber  were  rigidly  connected  to  a  piston/cylinder 
assembly  used  to  oscillate  the  fluid  inside  the  main  chamber.  The  platinum  wire  was  used 
as  the  heated  surface.  The  amplitude  and  frequency  of  the  oscillation  in  the  main  chamber 
was  varied  to  attempt  to  locate  an  optimum  value  with  a  small  amount  of  data.  A  local 
optimum  value  was  located  at  amplitude  of  0. 1 17mm  at  a  frequency  of  0.69F£z. 
C.  OBJECTIVES 

The  purpose  of  the  present  study  is  to  continue  the  work  begun  by  Turk  to 
investigate  the  effects  of  fluid  oscillation  on  the  onset  of  nucleate  boiling.  A  similar 
experimental  apparatus  with  an  enhanced  data  acquisition  system  was  used;  a  main 
chamber  filled  with  FC-72  housed  the  platinum  wire.  A  front  chamber  mounted  to  the 
main  chamber  and  serve  as  the  connection  with  the  piston  assembly.  One  major  change  is 
that  the  piston  assembly  is  not  rigidly  connected  to  the  front  or  main  chamber.   This  was 


changed  in  an  effort  to  eliminate  any  vibrations  due  to  the  piston  assembly  being 
transmitted  to  the  main  chamber  and  effecting  the  results.  An  attempt  will  be  made  to  find 
the  frequency  and  amplitude  of  oscillation  that  will  have  the  greatest  reduction  in  the 
excess  superheat.  Keeping  wall  superheat  as  low  as  possible  will  reduce  the  surface  heat 
flux  and  increase  the  reliability  of  the  electronic  equipment. 


H.  EXPERIMENTAL  APPARATUS 

A.         DESCRIPTION  OF  COMPONENTS 

The  experimental  apparatus  (Figure  1)  consists  of  four  main  parts:  the  main 
chamber  in  which  the  boiling  takes  place,  the  front  chamber,  the  adjustable  stroke  piston 
with  DC  motor  which  oscillates  the  fluid,  and  the  condenser  system  located  on  the  top  of 
the  main  chamber.  The  main  chamber  houses  the  platinum  wire  and  numerous 
thermocouples.  The  front  chamber  is  rigidly  mounted  to  the  main  chamber  and  serves  as  a 
connection  between  the  tubing  leading  to  the  piston/cylinder  oscillation  generator  and  the 
main  chamber.  Only  flexible  tubing  connects  the  front  chamber  and  the  adjustable  stroke 
piston.  The  condenser  system,  which  form  the  top  of  the  main  chamber  was  constructed 
for  this  experiment  and  is  designed  to  be  cooled  by  a  constant  temperature  water  bath. 

The  front  and  main  chambers  hold  the  dielectric  fluid  FC-72.  The  front  chamber 
regulates  the  motion  in  the  fluid,  which  is  caused  by  the  piston  and  transmits  the  motion  to 
the  main  chamber.  To  prevent  lose  of  FC-72  during  degassing,  all  seams  are  sealed  with 
Dow  Corning  3140  RTV  coating  and  Devcon  Corporation  Silicone  Adhesive  to  prevent 
leaks.  Also  at  the  clearance  between  the  piston  seal  on  the  front  chamber,  Teflon  tape  is 
utilized  to  prevent  leakage. 

1.  Main  Chamber  (Figure  2) 

The  main  chamber  is  a  Plexiglas  box  with  an  aluminum  heat  exchanger  serving  as 
the  top  cover.  The  inner  dimensions  of  the  main  chamber  are  152  x  60  x  152  mm.  (Figure 
2).    An  accurate  measurement  of  the  main  chamber  is  necessary  to  calculate  the  vertical 


displacement  of  the  fluid  in  the  main  chamber  as  a  function  volumetric  displacement  in  the 
piston/cylinder.  This  allows  a  direct  correlation  of  the  amplitude  of  the  induced  motion  of 
the  fluid  in  the  main  chamber  with  the  piston  stroke  length. 

A  Plexiglas  board,  which  holds  the  0.05mm  platinum  wire,  is  inserted  vertically 
down  the  center  of  the  main  chamber.  The  board  can  be  easily  placed  into  or  removed 
from  the  main  chamber  by  means  of  two  slots  on  the  sidewalls.  A  tube  at  the  back  wall  of 
the  main  chamber  is  used  to  house  the  necessary  wiring  for  the  platinum  wire  and  the 
numerous  thermocouples  installed  throughout  the  main  chamber.  All  of  the  wiring 
necessary  exits  via  this  tube  and  goes  either  to  the  connection  board  or  to  the  data 
acquisition  system. 

The  FC-72  needs  to  be  degassed  via  boiling  during  the  experiment.  The  FC-72  has 
is  capable  of  having  dissolved  air  up  to  48%  by  volume  in  solution.  [Ref.  3]  Dissolved  air 
in  the  FC-72  will  enhance  the  heat  transfer  properties  of  the  fluid  and  will  skew  the  results. 
Therefore,  bulk  heaters  have  been  installed  at  the  bottom  of  the  main  chamber. 
Thermocouples  have  been  installed  2cm  above  the  bottom  of  the  chamber  to  record  the 
fluid  bulk  temperature.  The  main  chamber  is  encased  in  insulation  to  reduce  the  amount 
of  heat  loss  to  the  atmosphere.  Another  important  feature  of  the  main  chamber  is  the  fluid 
entrance  port.  The  passage  from  the  front  chamber  to  the  main  chamber  is  an  oval 
opening  with  a  porous  plug  installed  to  evenly  smooth  the  effects  of  the  pump. 

The  main  chamber  consists  of  three  major  components:  the  bulk  heaters,  the 
condenser/heat  exchanger,  and  the  board  with  the  adjustable  spring  assembly. 
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a.  Heaters 

Three  strip  heaters  have  been  installed  at  the  bottom  of  the  main  chamber. 
The  heaters  are  used  to  provide  the  necessary,  bulk  heating  throughout  the  experiment. 
The  heaters  are  electric  powered  and  attached  to  a  thin  aluminum  plate  at  the  bottom  of 
the  main  chamber.  The  wiring  for  the  heaters  enters  the  chamber  via  the  tubing  at  the  rear 
of  the  unit.  During  the  experiment,  the  fluid  must  be  maintained  close  to  56°C,  which  is 
saturation  temperature  for  FC-72  at  atmospheric  pressure. 

The  dimensions  of  the  each  heater,  manufactured  by  MTNCO  Co.,  are  125 
x  1 1  mm.  The  heaters  are  in  parallel  connection  and  all  three  heaters  are  used  at  2.2  A  and 
20  V  during  the  degassing  process.  Two  heaters  are  used  during  the  experiment  at  0.6  A 
and  14  V.  For  that  reason  two  of  the  heaters  are  connected  to  the  power  supply  with 
toggle  switches.  The  heaters  that  remain  in  operation  throughout  the  experiment  are  the 
two  located  at  the  rear  of  the  chamber  and  behind  the  platinum  wire. 

b.  The  Board  and  Adjustable  Spring  Assembly  (Figure  3) 

The  0.05  mm  diameter  platinum  wire  and  the  spring,  which  is  used  to 
provide  tension  on  the  wire  at  all  times,  are  mounted  on  a  removable  Plexiglas  board  that 
is  inserted  into  the  main  chamber.  The  platinum  wire  has  been  soldered  to  the  brackets 
under  no  tension.  The  spring  assembly  can  be  adjusted  to  impose  the  proper  amount  of 
tension  for  the  conduct  of  the  experiment.  The  power  supply  wires  for  the  platinum  wire 
will  enter  and  exit  the  main  chamber  via  the  exit  tubing  in  the  rear  of  the  chamber. 


11 


c.  Condenser/Heat  Exchanger  (Figure  4) 

An  aluminum  heat  exchanger  serves  as  the  cover  for  the  main  chamber. 
This  heat  exchanger  condenses  the  FC-72  vapor  and  returns  it  to  the  chamber.  The  heat 
exchanger  is  secured  to  the  top  of  the  main  chamber  by  ten  screws  with  a  cork  gasket 
inserted  between  the  heat  exchanger  and  the  main  chamber.  The  condenser/heat 
exchanger  is  2.5cm  thick  plate  with  the  bottom  0.6cm  solid  aluminum.  Six  grooves, 
1.8cm  deep,  were  carved  into  the  plate  to  serve  as  passages  for  the  cooling  liquid  to  flow. 
A  Plexiglas  cover  is  mounted  on  the  top  of  the  heat  exchanger  with  a  cork  gasket  inserted 
between  the  cover  and  the  heat  exchanger.  The  Plexiglas  cover  contains  the  cooling 
liquid's  enter  and  exit  fitting.  A  pump,  which  is  part  of  the  constant  temperature  bath, 
circulates  the  cooling  liquid,  which  is  distilled  water.  The  constant  temperature  bath, 
(RTE-111  Refrigerated  Bath/Circulator,  Model  No.  134103201100)  is  manufactured  by 
the  NESLAB  Instruments  Inc.  It  has  a  temperature  range  of  -25  C  to  150  C  with  a 
temperature  stability  of +/-  0.01  C.  This  unit  has  a  1.9  gallon  reservoir  and  the  cooling 
liquid  is  transported  in  0.25  inch  inner  diameter  tygon  tubing.  The  bath  has  a 
microcomputer  temperature  control  system  that  will  continually  adjust  the  system  to 
maintain  pump  outlet  temperature  at  the  set  temperature.  Thermocouples  are  installed  at 
the  inlet  and  outlet  fitting  to  the  heat  exchanger  to  monitor  cooling  water  temperature 
change  across  the  heat  exchanger. 

2.  Front  Chamber  (Figure  5) 

The  front  chamber  has  been  added  to  the  main  chamber  to  assist  in  removing  any 
irregularities  from  the  flow  prior  to  entry  in  the  main  chamber.  A  porous  plug  is  installed 
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over  the  opening  into  the  main  chamber  to  serve  as  a  flow  dampener.  These  irregularities 
could  affect  the  oscillation  of  the  fluid  inside  the  main  chamber.  Any  interference  with  the 
amplitude  of  oscillation  will  alter  the  results  of  the  study. 

3.  DC  Motor  and  Piston/Cylinder  Assembly  (Figure  6) 

For  the  experiment,  a  piston/cylinder  assembly  has  been  fabricated  to  act  as  the 
fluid  oscillation  generator.  The  piston  has  an  outer  diameter  of  9.5mm  (0.375  inches)  and 
has  been  designed  to  create  an  oscillatory  motion  in  the  fluid  to  assist  in  the  removal  of 
bubbles  from  the  platinum  wire. 

The  piston/cylinder  assembly  is  not  rigidly  connected  to  the  main  chamber.  The 
fluid  oscillation  is  transported  through  2.75  feet  of  flexible  rubber  tubing  that  has  been 
coiled  to  assist  in  isolating  the  main  chamber.  This  has  been  implemented  to  attempt  to 
prevent  any  vibrations  from  the  DC  motor  affecting  the  main  chamber's  fluid  motion. 

The  frequency  and  the  amplitude  of  the  oscillatory  motion  in  the  main  chamber  are 
determined  from  the  variable  stroke  and  frequency  of  the  piston.  A  slot  on  the  rotating 
disk  connected  to  the  DC  motor  achieves  the  variable  piston  stroke.  By  adjusting  the 
placement  of  the  connection  point  in  the  slot,  the  stroke  length  can  be  adjusted  as 
required.  The  stroke  length  will  be  measured  using  a  depth  micrometer  with  a  reference 
point  established  on  the  piston  assembly  foundation.  By  varying  the  stroke  length,  the 
amount  of  FC-72  moving  from  the  piston/cylinder  assembly  to  the  main  chamber  is  varied. 
The  stroke  length  of  the  piston/cylinder  assembly  controls  the  amplitude  of  oscillation 
inside  the  main  chamber.  A  27  volt,  1  Amp  DC  motor,  achieves  the  variable  frequency. 
The  rotational  speed  is  measured  by  a  micro  switch  tachometer. 
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B.  INSTRUMENTATION  (Figure  7) 

The  main  goal  of  this  experiment  is  to  observe  changes  in  the  boiling  curve  of  FC- 
72  by  obtaining  the  surface  temperature  of  the  platinum  wire  at  various  surface  heat  flux 
values.  For  that  reason  the  platinum  wire  temperature  vs.  resistance  dependence  is 
calibrated.  The  results  of  the  platinum  wire  calibration  is  a  linear  temperature  vs. 
resistance  relationship.  The  platinum  wire  is  connected  in  series  to  a  precision  resistor  of 
known  resistance  in  series.  With  this  arrangement,  measuring  the  voltage  drop  across  the 
precision  resistor  allows  the  current  in  the  platinum  wire  to  be  calculated.  Next  the 
voltage  drop  across  the  platinum  wire  is  measured  and  the  resistance  of  the  wire  is  then 
easily  calculated.  With  the  known  resistance  value,  the  surface  temperature  of  the 
platinum  wire  can  be  obtained  from  the  calibration.  Since  the  voltage  drop  and  current  in 
the  platinum  wire  are  known,  their  product  yields  the  heat  transfer  rate  and  division  by  the 
surface  area  of  the  platinum  wire  will  yield  the  surface  heat  flux.  [Ref  9] 

The  bulk  temperature  of  the  fluid  and  the  condenser  temperature  are  measured  by 
thermocouples  inside  the  main  chamber.  Some  additional  instrumentation  is  required  to 
obtain  the  frequency  of  the  DC  motor  and  the  piston.  Figure  6  is  the  wiring  diagram  of 
power  supplies  and  data  acquisition  system  without  thermocouples. 

1.  Platinum  Wire 

The  platinum  wire  is  a  product  of  Goodfellow  Advanced  Materials  Ltd.  of 
Cambridge,  England.  It  has  99.9%  purity  and  a  nominal  diameter  of  0.05mm.  The  active 
length  of  the  platinum  wire  is  102.8mm.  Observations  under  a  scanning  electron 
microscope  of  a  piece  of  wire  from  the  same  spool  revealed  a  diameter  of  51.lu.rn.  When 
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mounting  the  wire,  very  careful  handling  of  the  wire  is  required  due  to  its  small  diameter, 
especially  when  moving  the  wire  into  the  main  chamber  and  when  removing  the  wire  from 
the  calibration  bath. 

2.  Power  Supplies 

In  this  set  of  experiments,  three  DC  power  supplies  are  used;  one  for  the  platinum 
wire,  one  for  the  bulk  heaters,  and  one  for  the  piston  assembly  DC  motor.  These  power 
supplies  are: 

HP  6289A  DC  Power  Supply  0-50V,  0-1. 8A  with  HP  59501B  Isolated 
DAC/Power  Supply  Programmer.  (HP  6289A  provided  voltage  to  the  platinum  wire  and 
the  voltage  value  is  controlled  by  the  Power  Supply  Programmer  and  the  computer 
program.) 

HP  6268 A  DC  Power  Supply  0-24 V,  0-1 2 A  (to  power  the  bulk  heaters.) 

HP  6289 A  DC  Power  Supply  0-50 V,  0-1. 8 A  (to  power  the  DC  motor.) 

3.  Data  Acquisition  System 

All  the  temperature  and  voltage  measurements  are  gathered  by  HP  3 852 A  Data 
Acquisition  and  Control  Unit  via  the  Unitek  desktop  computer  with  Pentium  processor. 
The  HP  3 8 52 A  Data  Acquisition  system  is  able  to  convert  temperature  data  directly  from 
the  thermocouple  to  degrees  Celsius.  The  cards  that  are  used  in  the  data  acquisition 
system  are: 

1.  Thermocouple  card  44708  A  20  channel  T-couple  Acquisition 

2.  Voltage  card  44705A  20  channel  Guarded  Acquisition 
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4.  Frequency  Measurements 

A  5001  Universal  Counter-Time  will  be  used  to  measure  the  period.     The 
frequency  will  be  obtained  from  this  measurement. 
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Figure  1.  Diagram  of  Experimental  System. 
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Figure  2.  Main  Chamber. 
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Figure  3.  Board  and  Adjustable  Spring  Assembly. 
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Figure  4.  Condenser/Heat  Exchanger. 
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Figure  5. 


Front  Chamber. 
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Figure  6.  DC  Motor  and  Piston/Cylinder  Assembly. 
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Instrumentation  Schematic  Diagram. 
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m.  EXPERIMENTAL  PROCEDURE 

A.  SYSTEM  PREPARATIONS 

Prior  to  each  experimental  run,  the  following  preparatory  steps  must  be  taken: 

1 .  Check  the  fluid  level  of  the  main  chamber  and  the  tube.  The  fluid  level  in  the 
main  chamber  needs  to  be  higher  than  the  platinum  wire  by,  and  no  higher  than,  2.5cm 
above  the  top  of  the  platinum  wire's  Plexiglas  mount.  The  fluid  level  in  the  tube  should  be 
no  higher  than  the  fluid  level  in  the  main  chamber.  If  necessary,  add  more  FC-72  via  the 
tube  in  the  rear  of  the  main  chamber. 

2.  Adjust  the  piston  stroke  length  to  the  desired  setting.  Ensure  the  adjusting  bolt 
is  securely  fastened. 

3.  Turn  on  the  HP3852A  Data  Acquisition  System  and  the  Unitek  desktop 
computer. 

4.  Ensure  the  heat  exchanger  is  securely  fastened  to  the  top  of  the  main  chamber 
and  the  constant  temperature  bath  reservoir  level  is  at  the  high  mark. 

5.  The  main  chamber  must  be  fully  insulated  during  degassing  and  during  the 
experiment. 

B.  DEGASSING 

Once  system  preparations  are  completed,  the  following  steps  are  to  be 
accomplished  to  ensure  proper  degassing  of  the  FC-72: 

1.  Power  all  three  bulk  heaters  using  the  HP  6268 A  DC  Power  Supply  at  20  V 
and  2.2  A.  Be  sure  that  the  top  two  toggle  switches,  which  are  on  the  connection  board, 
are  on. 

25 


2.  Turn  on  the  constant  temperature  bath.  Set  the  bath  temperature  at  20  °C. 
Check  system  fully  to  ensure  that  no  leaks  are  present. 

Watch  the  fluid  temperature  increase  and.  the  temperature  drop  across  the  heat 
exchanger  steady  out  at  +/-  0.3  °C.  The  fluid  temperature  should  rise  up  to  the  saturation 
temperature  of  56  °C.  At  this  temperature  degassing  will  commence  with  pool  boiling. 
The  complete  degassing  process  takes  approximately  3  hours  to  complete. 
C.        EXPERIMENTAL  PROCEDURE 

The  following  steps  are  to  be  followed  before  commencing  the  experiment: 

1 .  Turn  off  the  bulk  heater  using  the  toggle  switch  at  the  top  of  the  connection 
board. 

2.  Decrease  the  power  of  the  remaining  heaters  to  12  V  and  0.6  A.  Wait  until  the 
temperature  drops  down  to  a  thermal  equilibrium  point. 

3.  Turn  on  the  5001  Universal  Counter-Timer. 

4.  Power  DC  motor  by  turning  on  the  HP  6289A  DC  Power  Supply.  Adjust  the 
DC  motor  voltage  to  obtain  the  desirable  frequency.  Verify  the  desired  frequency  has 
been  obtained  via  the  Universal  Counter-Timer  display.  Once  properly  set;  use  toggle 
switch  on  connection  board  to  stop  motor. 

5.  Turn  on  the  HP6289A  DC  Power  Supply  and  HP59501B  Isolated  DAC/Power 
Supply  to  power  the  platinum  wire  and  the  precision  resistor. 

6.  Load  the  QBASIC  computer  program  named  Tuitel.exe 

7.  Run  first  four  power  setting  without  running  the  DC  motor  for  the  piston. 

8.  Save  each  data  set  to  a  3.5  inch  disk  inserted  in  drive  A. 
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9.    At  the  end  of  the  fourth  power  setting,  turn  on  the  DC  motor  for  the  piston  at 
the  connection  board  via  the  toggle  switch. 
D.  DATA  ACQUISITION  PROGRAM 

Data  acquisition  is  accomplished  using  a  QBASIC  program  called  Tuitel.exe. 

1 .  The  data  acquisition  program  executes  these  operations: 

Sets  up  communication  with  the  HP3852A  and  HP6289A  DC  Power 
Supply. 

Controls  power  settings  incrementally  by  the  help  of  HP59501B  Power 
Supply  Programmer. 

Reads  thermocouple  temperatures  after  the  HP3852A  Data  Acquisition 
System  converts  the  thermocouple  voltages  to  temperature  values  (°C). 

Reads  voltage  drops  across  the  precision  resistor  and  the  platinum  wire. 

Stores  all  recorded  values  into  a  data  file  either  to  the  hard  drive  or  a  3.5 
inch  diskette. 

2.  For  each  power  setting  these  measurements  will  be  taken  once  a  steady  state 
has  been  reached. 

Bulk  temperature  (four  thermocouples). 

Ambient  temperature  (two  thermocouples). 

Heat  exchanger  inlet  temperature  (one  thermocouple). 

Heat  exchanger  outlet  temperature  (one  thermocouple). 

Platinum  wire  voltage  drop. 

Precision  resistor  voltage  drop. 
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3.  The  first  five  steps  of  power  setting;  the  voltage  increment  is  0.5  V.  After  step 
five,  the  increment  will  be  0.4  V  for  the  next  four  power  settings.  After  step  eight  the 
increment  will  be  0.2  V  until  step  22.  At  the  end  of  the  22nd  reading,  the  power  setting 
will  start  to  decrease  by  0.4  V  per  step.  After  the  end  of  the  31st  reading,  the  power 
setting  will  be  decreased  by  0.5  V  per  step  and  will  end  after  the  35th  reading  at  1 .0  V. 

4.  Repeat  the  same  steps  for  all  other  piston  stroke  and  frequency  configurations. 
E.  DATA  REDUCTION 

The  following  calculations  are  used  to  reduce  the  raw  data  obtained  by 
experiments. 

1 .    Condenser  temperature: 


T      m  =  (T0+T}+T2+T3) 

A  cond  V   ^-v 


To  through  T3  are  thermocouple  temperatures. 

2.  Fluid  bulk  temperatures: 

T     m_  (T4+T5+T6+T7) 

4 

T4  through  T4  are  thermocouple  temperatures. 

3.  Current  of  platinum  wire  and  precision  resistor: 


IptWire(A)  -  I2Q  - 


_  r  2Q 


^2Q 


Where  R2Q  =  2.03Q  (Precision  Resistor) 
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4.  Platinum  wire  resistance  value: 


RptWire(£*)  — 


PtWire 
*  PtWire 


5.  Platinum  wire  surface  temperature: 

Tsurf  (°C)  —  8  *  Rpt.Wire  +  Ty=o 

where: 

s  =  49.689  (slope  of  calibration  curve) 

Ty=0  =  -253.836  (y  axis  intercept  of  calibration) 

6.  Fluid  film  temperature: 

7L,+7L 


Tfilm(°C)  - 


'  surf  bulk 


2.0 
7.    Heat  flux  from  platinum  wire: 


I         *V 
q"(W/m  )  = 


A 

■^PtWire 


8.    Piston  displacement: 


Apiston(mm3)  -  —  *  D2Plston  *  LSlroke 
4 

9.  Amplitude  of  oscillation  in  the  main  chamber: 

AAmplitude(mm)  =  '■L-fL 

a*b 
Where  a  and  b  are  the  inner  length  and  width  of  the  main  chamber. 
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9.  Frequency  calculation: 

1 


/Hz) 


*  Penod  \S) 


A  MATLAB  program  was  written  to  evaluate  the  gathered  data.  The  program 
was  written  in  three  parts.  The  first  part  would  retrieve  the  gathered  data  from  the 
appropriate  data  file  and  apply  the  temperature  corrections  for  the  thermocouples.  The 
second  part  would  calculate  the  heat  flux  based  on  the  gathered  data  and  the  calculated 
surface  temperature  of  the  platinum  wire.  Next  the  program  would  calculate  the  natural 
convection  using  the  recommended  correlation  by  Kuehn  and  Goldstein  [Ref  2]  for 
natural  convection  around  small  horizontal  cylinders.  Then  both  sets  of  results  are  plotted 
versus  the  temperature  superheat  value.  For  fluid  properties,  Egger's  thesis  [Ref.  11]  is 
used. 
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IV.        RESULTS  AND  DISCUSSION 

In  the  present  study,  the  boiling  curve  of  FC-72  is  determined  in  a  non-oscillating 
bulk  fluid.  Changes  in  the  boiling  curve  of  FC-72,  specifically  the  value  of  the  wall 
superheat  at  the  initiation  of  nucleate  boiling,  are  observed  when  the  bulk  fluid  is 
oscillating  at  different  amplitudes  and  frequencies.  All  experiments  are  conducted  at 
atmospheric  pressure,  21°C  ambient  temperature,  and  the  same  platinum  wire  of 
0.0511mm  diameter  is  used. 

The  results  of  these  experiments  are  presented  in  the  plots  of  wire  heat  flux  vs. 
wall  superheat  (Tsurf-  T^t).  Figures  8  through  12  are  a  representative  sample  of  the  data 
gathered  during  the  experiments.  The  remaining  graphs  are  included  in  APPENDIX  E. 
Figure  8  is  the  boiling  curve  of  FC-72  without  any  oscillation  in  the  bulk  fluid.  The 
hysteresis  between  the  increasing  and  the  decreasing  heat  flux  shows  the  amount  of  wall 
superheat  necessary  for  the  initiation  of  nucleate  boiling.  The  surface  temperature  of  the 
wire  decreased  abruptly  with  the  initiation  of  nucleate  boiling.  Figure  10  is  the  boiling 
curve  of  FC-72  with  an  amplitude  of  oscillation  of  0.105mm  and  an  oscillation  frequency 
of  1.4Hz.  The  decrease  in  the  magnitude  of  the  hysteresis  is  obvious  when  compared  to 
Figure  8. 
A.         BOILING  CURVE  OF  FC-72  (BOILING  WITHOUT  OSCILLATION) 

To  study  the  effects  of  oscillation  on  the  boiling  curve  of  FC-72,  it  is  necessary  to 
establish  the  boiling  curve  of  FC-72  without  oscillation  of  the  bulk  fluid.     The  boiling 
curve  of  FC-72  under  non-oscillating  condition  will  be  the  base  line  for  comparison  for 
subsequent  experiments.    Three  experiments,  Experiments  4,  6,  and  25,  were  conducted 
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under  non-oscillating  conditions.  The  results  of  Experiment  6  are  shown  in  Figure  8. 
Experiments  4  and  25  results  are  in  agreement  with  Experiment  6  and  are  included  in 
APPENDIX  E.  The  maximum  wall  superheat  (Tsulf  -  Tsat)  value  prior  to  nucleate  boiling 
was  45 °C  for  this  part  of  the  experiment.  Maximum  wall  superheat  values  for  each 
experiment  in  this  section  are  listed  in  Table  1 .  The  fluid  bulk  temperature  was  maintained 
between  51°C  and  52.5°C  for  this  experiment.  In  the  natural  convection  phase  of  the 
boiling  curve  of  FC-72,  the  data  is  in  agreement,  within  uncertainty  limits,  with  the 
correlation  suggested  by  Keuhn  and  Goldstein  [Ref.  3]  for  natural  convection  heat  transfer 
around  small  diameter  cylinders.  The  platinum  wire  surface  temperature  just  after  the 
initiation  of  nucleate  boiling  is  19.9%  less  than  the  maximum  wall  superheat  value.  The 
percent  drop  in  platinum  wire  surface  temperature  just  after  the  initiation  of  nucleate 
boiling  is  in  agreement  with  Egger's  results.  [Ref.  11]. 


Experiment 
Number 

Wall  superheat  necessary  to 
initiate  nucleate  boiling 
':J     <°C) 

4 

45.3 

6 

45.0 

25 

45.7 

Table  1. 


The  Comparison  of  Maximum  Wall  Superheat 
Values  for  Three  Non-Oscillating  Experiments. 


B. 


OSCILLATION  EFFECTS  ON  THE  BOILING  CURVE 


Four  different  oscillation  amplitudes  in  the  fluid  surrounding  the  platinum  wire 
were  employed  during  the  study.  Turk  [Ref.  10]  reported  a  local  optimum  oscillation 
amplitude  of  0.117mm,  which  was  used  in  this  study  as  a  starting  point.  Oscillation 
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amplitudes   of  0.105,   0.125,   and  0.252mm  respectively  were   also  used   at  various 

frequencies  ranging  from  0.3  to  1.5Hz  respectively. 

1.  Amplitude  and  Frequency  Effects  on  the  Boiling  Curve 

Table  1  represents  the  maximum  wall  superheat  values  observed  prior  to  the 

initiation  of  nucleate  boiling  for  the  non-oscillating  and  oscillating  FC-72. 

a.  A  mplitude  of  0.117  mm 

The  oscillation  amplitude  of  0.117mm  was  identified  as  a  local  optimum 
value  for  oscillation  amplitude  in  Turk's  report.  [Ref.  10].  The  experiment  was  conducted 
using  oscillation  frequencies  of  0.69,  1.0,  1.2,  and  1.4Hz  respectively.  The  lowest  wall 
superheat  value  to  initiate  nucleate  boiling  was  29.1°C,  which  is  a  33%  drop  from  the 
boiling  curve  of  FC-72,  and  was  obtained  at  an  oscillation  frequency  of  0.69Hz.  The  heat 
flux  vs.  wall  superheat  plot  of  this  frequency  is  Figure  9.  Wall  superheat  values  at  various 
frequencies  are  listed  in  Table  2.  The  best  results  in  the  wall  superheat  drops  are  achieved 
at  this  amplitude.  As  the  frequency  of  oscillation  was  increased  the  wall  superheat 
necessary  to  initiate  nucleate  boiling  also  increased  with  a  maximum  value  of  37.2°C  at  1.4 
Hz.  Even  the  maximum  value  of  wall  superheat  is  a  20.9%  drop  in  wall  superheat 
compared  to  the  non-oscillating  values.  Wire  surface  temperature  and  heat  flux  at  various 
frequencies  are  listed  in  Table  2.  The  remaining  heat  flux  vs.  wall  superheat  plots  are 
included  in  APPENDIX  E. 

b.  Amplitude  of  0.105  mm 

The  oscillation  amplitude  of  0.105mm,  which  is  approximately  two  wire 
diameters  in  size,  produced  a  decrease  in    the    wall  superheat  value  necessary  to  initiate 
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Amplitude 

0.105  mm 

0.117mm 

0.125  mm 

Frequency 
(Hz) 

Wire 
Temp 

CO 

Bulk 
Temp 

(°Q 

Heal 

Ffux 

(W/nT) 

Wire 

Temp 

CO 

Bulk 
Temp 

CO 

Heat 
Flux 

(W/m2) 

Wire; 

Temp; 

CO 

Bulk 

Temp 

CO 

Heat 
Flux 

(W/m2) 

No  Stroke 

101.1 

52.9 

81,300 

101.1 

52.9 

81,300 

101.1 

52.9 

81,300 

0.6 

93.3 

53.3 

84,072 

— 

~ 

— 

98.5 

52.3 

90,295 

0.69 

- 

— 

— 

85.2 

52.2 

76,082 

— 

— 

— 

1.0 

90.1 

53.1 

82,947 

86.0 

51.3 

92,474 

91.5 

51.8 

75,443 

1.2 

87.9 

52.7 

68,104 

90.0 

52.4 

75,836 

91.7 

52.1 

75,895 

1.4 

86.2 

52.7 

100,000 

92.5 

53.6 

75,313 

98.7 

53.2 

90,309 

Table  2.  Superheat  Values  for  Oscillating  FC-72.  For  No-Stroke 

Case  the  Results  of  Experiment  6  is  Used. 

nucleate  boiling  as  the  oscillation  frequency  was  increased.    The  oscillation  frequency  of 

1.4Hz  produced  a  wall  superheat  value  of  30.5°C.    The  heat  flux  vs.  wall  superheat  plot 

for  this  frequency  is  Figure  10.    This  is  an  improvement  from  a  value  of  37.2°C  for  wall 

superheat  that  was   achieved  at  an  oscillation  frequency  of  0.6F£z.      Wire   surface 

temperature  and  heat  flux  for  various  frequencies  are  listed  in  Table  2.  The  remaining  heat 

flux  vs.  wall  superheat  plots  for  this  amplitude  are  included  in  APPENDIX  E. 

c.  Amplitude  of  0.125  mm 

The  oscillation  amplitude  of  0.125mm,  which  is  approximately  2.5  wire 
diameters,  produced  an  erratic  decrease  in  wall  superheat  values  necessary  to  initiate 
nucleate  boiling.  A  minimum  wall  superheat  value  of  35.4°C  was  achieved  at  an  oscillation 
frequency  of  1.0  and  35.8°C  at  1.2Hz  respectively.  A  higher  wall  superheat  value  was 
recorded  at  oscillation  amplitudes  of  0.6Hz  and  1.4Hz.  The  heat  flux  vs.  wall  superheat 
plot  for  a  frequency  of  1.4  Hz  is  Figure  11.  At  an  amplitude  of  0.125mm,  the  maximum 
value  of  wall  superheat  is  an  8%  drop  in  wall  superheat  compared  to  the  non-oscillating 
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value.    Wall  superheat  values  at  various  frequencies  are  listed  in  Table  2.    The  remaining 
heat  flux  vs.  wall  superheat  plots  for  this  amplitude  are  included  in  APPENDIX  E. 
C.         AMPLITUDES  INDEPENDENT  OF  FREQUENCY 

At  an  amplitude  of  oscillation  of  0.252mm,  which  is  approximately  5  wire 
diameters,  the  wall  superheat  value  necessary  to  initiate  nucleate  boiling  appears  to  be 
independent  of  oscillation  frequency.  Oscillation  frequencies  of  0.3,  1.0,  and  1.5  Hz 
respectively  were  employed  and  all  three  produced  a  wall  superheat  value  of  30.5°C 
±0.5°C.  The  heat  flux  vs.  wall  superheat  plot  for  a  frequency  of  1.0Hz  is  Figure  12.  Wall 
superheat  values  at  various  frequencies  are  included  in  Table  3.  The  remaining  heat  flux 
vs.  wall  superheat  plots  for  this  amplitude  are  included  in  APPENDIX  E. 


Amplitude 

0.252  mm 

Frequency 
<Hz)       :  ; 

Wire 
Temp 

m 

Bulk 

Temp 

<°C) 

Heat 

Flux 

(W/m2) 

0.3 

87.3 

51.1 

69,542 

1.0 

86.95 

51.5 

69,874 

1.5 

86.6 

51.1 

83,334 

D. 


Table  3 .  Superheat  Values  for  Oscillating  FC-72  at  an 

Amplitude  of  0.252  mm. 

FURTHER  DISCUSSIONS 


A  total  of  32  experiments  have  been  conducted  to  achieve  the  boiling  curve  of  FC- 
72  at  various  fluid  oscillation  amplitudes  and  frequencies.  The  effects  of  the  fluid 
oscillation  amplitude  and  frequencies  on  the  wall  superheat,  which  is  necessary  to  initiate 
nucleate  boiling,  have  been  observed.  Figure  13  shows  the  maximum  wall  superheat 
values  required  to  initiate  nucleate  boiling  at  four  different  amplitudes  and  various 

frequencies.   An  overall  result  of  the  oscillatory  condition  is  a  certain  decrease  in  required 
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wall  superheat  for  boiling  incipience.  The  oscillation  amplitude  of  0.117mm  at  an 
oscillation  frequency  of  0.69Hz  is  the  most  effective  combination  attempted  during  this 
study. 

The  reduction  of  wall  superheat  necessary  to  initiate  nucleate  boiling  can  be 
attributed  to  the  way  bubbles  forming  at  the  nucleation  sites  depart  from  the  heated 
surface.  When  the  bubble  outgrows  the  nucleation  site,  it  starts  to  grow  bigger  on  the 
wire  surface.  In  highly  wetting  fluids  like  FC-72,  the  bubbles  can  not  depart  the  surface 
easily  due  to  low  surface  tension  of  the  fluid  so  the  surface  temperature  of  the  heated 
surface  continues  to  increase.  The  bubble  departing  process  from  nucleation  sites  and  the 
heated  surface  is  accelerated  due  to  the  oscillation.  The  oscillations  assist  the  bubbles 
overcoming  the  surface  tension  effects  of  highly  wetting  FC-72  by  mechanically  removing 
the  bubbles.  As  was  expected,  some  amplitudes  and  frequencies  of  oscillation  were  more 
effective  than  others  at  reducing  wall  superheat  values  necessary  to  initiate  nucleate 
boiling.  It  is  very  important  to  consider  that  in  every  case  involving  fluid  oscillation  the 
value  of  required  wall  superheat  was  lower  than  in  the  non-oscillating  fluid  case. 
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In  Figures  8  to  32  These  Symbols  Are  Used: 
X  Increasing  Heat  Flux  (TSUIf  -  Tsat) 

O  Decreasing  Heat  Flux  (Tsurf  -  T^t) 

+  Increasing  Heat  Flux  (Tbuik) 

Decreasing  Heat  Flux  (Tbuik) 

Kuehn-Goldstein  Correlation  (Natural  Convection) 
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Heat  flux  vs.  Wall  superheat 
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Figure  8.  Boiling  Curve  of  FC-72. 
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Heat  flux  vs.  Wall  superheat 


Exp  No  14    Freq  =  0.69Hz    Amp  =  0.1 17  mm 


li.  10" 

O) 

X 


10 


10 


20 


30  40  50  60 

Wall  Superheat  deg  C 


Figure  9.  Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude  0. 1 17mm 

And  a  Frequency  of  0.69Hz. 


39 


10 


Heat  flux  vs.  Wall  superheat 


Exp  No.  22    Freq  =  1 .4  Hz  Amp  =  0.1 05  mm 
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Figure  10.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude  0. 105mm 
And  a  Frequency  of  1  4Hz. 
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Heat  flux  vs.  Wall  superheat 
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Figure  1 1 .        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude  0. 125mm 
And  a  Frequency  of  1 .4Hz. 
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Heat  flux  vs.  Wall  superheat 
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Figure  12.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude  0.252mm 
And  a  Frequency  of  1  5Hz. 
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Frequency  vs  Wall  Superheat 


Frequency  (Hz) 


Figure  13.        Results  of  Experiments  -  Heat  Flux  vs.  Wall  Superheat  at 
Various  Amplitude  and  Frequencies  of  Oscillation. 
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V.         CONCLUSIONS  AND  RECOMMENDATIONS 

A.         CONCLUSIONS 

The  present  study  investigated  boiling  of  degassed  FC-72  in  an  oscillating  fluid 
condition.  An  enhancement  of  the  transition  to  nucleate  boiling  has  been  confirmed  by  the 
results  achieved  in  this  study.  All  attempted  oscillation  amplitudes  and  frequencies 
changed  the  boiling  curve  of  highly  wetting  FC-72.  The  largest  required  wall  superheat 
value  (Tsurf  -  Tsat)  necessary  to  initiate  nucleate  boiling  was  recorded  in  the  non-oscillating 
fluid  experiment.  At  some  amplitudes  the  required  wall  superheat  decreased  with 
increasing  oscillation  frequency  while  other  cases  had  the  required  wall  superheat  increase 
with  increasing  oscillation  frequency.  The  amplitude  of  0.117mm  and  the  frequency  of 
0.69Hz  provided  the  lowest  wall  superheat  temperature  for  nucleate  boiling  at  29.1°C. 
The  oscillation  amplitude  of  0.252mm  appeared  to  be  independent  of  oscillation  frequency 
since  a  required  wall  superheat  value  of  30.5°C  ±0.5°C  was  achieved  at  three  different 
frequencies. 

Even  though  the  experiments,  which  are  conducted  in  this  study,  proved  to  be 
reproducible  with  uncertainty  limits,  in  the  nucleate  boiling  phase  the  wall  superheat  value 
displayed  an  erratic  behavior.  This  was  expected  due  to  the  unsteady  nature  of  the  boiling 
phenomena.  No  two  experiments  will  be  able  to  exactly  reproduce  the  same  results  due  to 
the  haphazard  nature  of  boiling  at  the  bubble  level. 

Here  it  is  possible  to  conclude  that  the  required  wall  superheat  value  necessary  to 

initiate  nucleate  boiling  can  be  decreased  with  oscillation  in  the  bulk  fluid.    The  reason  can 

be  attributed  to  the  oscillating  fluid  removing  the  bubbles  from  the  heated  surface  prior  to 
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their  normal  departure  time.  This  makes  nucleation  sites  available  for  future  bubble 
growth  at  a  quicker  rate  than  normal.  These  nucleation  sites  are  now  acting  in  a  more 
efficient  manner  and  are  able  to  activate  other  sites  as  well.  This  will  enhance  the 
transition  to  nucleate  boiling  and  will  reduce  the  required  wall  superheat  value.  These 
results  are  in  agreement  with  the  results  in  Turk  [Ref.  10]  and  Egger  [Ref.  11]. 
B.         RECOMMENDATIONS 

In  this  study  only  four  amplitudes  of  oscillation  were  tested  at  selected  frequencies. 
For  future  study  a  more  detailed  frequency  and  amplitude  map  should  be  created.  Also,  at 
the  oscillation  amplitude  of  0.252mm,  the  required  wall  superheat  necessary  to  initiate 
nucleate  boiling  appeared  to  be  independent  of  frequency.  Continue  experiments  at 
oscillation  amplitudes  between  the  last  two  amplitudes  tested,  0.125  and  0.252mm 
respectively,  to  attempt  to  locate,  if  one  exists,  an  oscillation  amplitude  where  wall 
superheat  effects  are  truly  independent  o  frequency. 

During  the  study  numerous  experimental  runs  needed  to  be  abandoned  due  to 
leaks  in  the  main  chamber.  The  highly  wetting  FC-72  penetrated  all  seams  in  the 
structure,  including  seams  coated  with  sealant.  A  new  experimental  apparatus  needs  to  be 
fabricated  using  a  one  piece  mold  for  the  main  chamber.  Since  the  piston/cylinder 
assembly  is  no  longer  rigidly  attached  to  the  main  chamber,  the  front  chamber  can  be 
smaller  and  be  fabricated  as  part  of  the  main  chamber. 
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APPENDIX  A.  CALIBRATION  OF  PLATINUM  WIRE 

For  this  study  a  platinum  wire  is  used  as  a  temperature  reading  device.  The 
purpose  of  this  calibration  exercise  is  to  establish  a  relationship  between  the  resistance  of 
the  platinum  wire  and  its  surface  temperature.  This  relationship  is  known  to  be  linear. 
The  amount  of  change  in  the  resistance  of  the  platinum  wire  can  be  calibrated  to  measure 
temperature.  Turk  [Ref.  10]  used  this  same  process. 

Nikuyama  [Ref.  11]  used  the  platinum  wire  as  a  temperature  reading  device  in  the 
experiments  to  determine  the  boiling  regime  of  water.  He  used  a  platinum  wire  to 
identify  different  regimes  of  pool  boiling  of  water.  The  temperature  of  the  wire  was 
determined  from  the  knowledge  of  the  manner  in  which  its  electrical  resistance  varied 
with  temperature.  [Ref.  10]. 
A.         SEM  ANALYSIS  OF  PLATINUM  WIRE 

The  platinum  wire,  from  the  same  sample  which  is  used  during  the  experiments, 
is  analyzed  in  Scanning  Electron  Microscope  under  two  magnifications.  The 
magnification  of  666x,  as  it  can  be  seen  in  Figure  14a,  yields  that  the  diameter  of  the 
platinum  wire  in  0.05113mm.  The  magnification  of  1270x  (Figure  14b)  shows  the 
lengthwise  line  formations  created  during  the  manufacturing  process.  The  surface  of  the 
wire  is  relatively  smooth  in  appearance  with  the  exception  of  one  raised  ridge  line 
extending  the  length  of  the  sample.  Several  small  surface  defects  of  approximately  1-2 
(im  are  diameters  are  present  on  the  surface  of  the  wire.  The  defects  are  dispersed 
throughout  the  surface  of  the  wire.  These  are  possible  nucleation  sites  for  nucleate 
boiling. 
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Figure  14a.       SEM  Picture  of  Sample  Platinum  Wire  (666x). 


Figure  14b.    SEM  Picture  of  Sample  Platinum  Wire  (1270x). 
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B.         CALIBRATION  PROCEDURE  IN  THE  CALIBRATION  BATH 

During  the  calibration  in  the  bath  the  following  equipment  was  used: 

1.  HP  3852A  Data  Acquisition  System  with  installed  Ohmmeter. 

2.  Rosemont  Engineering  Company  Model  162C  Serial  No.  985  Platinum 
Resistance  Thermometer  (PRT)  (to  measure  the  calibration  bath  temperature. 
According  to  Reference  11,  PRT  has  been  calibrated  in  comparison  with  a 
platinum  resistance  temperature  standard,  which  is  done  by  the  National 
Bureau  of  Standards.) 

3.  Rosemont  Engineering  Company  Model  923B  Power  Supply  (to  control  bath 
temperature  with  Neslab  Endocal  Refrigerated  circulation  bath) 

4.  Rosemont  Engineering  Company  Model  920A  Commutating  Bridge  (to 
measure  PRT  resistance  for  temperature  measurement) 

5.  Rosemont  Engineering  Company  Model  913  A  Calibration  Bath  (  fluid  in  the 
bath  is  ethylene  glycol) 

The  0.05  mm  platinum  wire  has  been  mounted  on  a  Plexiglas  board  and  the  spring 
assembly  applies  tension  in  order  to  keep  the  wire  straight  throughout  the  calibration 
process.  The  board  and  spring  assembly  is  immersed  into  the  Rosemont  Engineering 
constant  temperature  bath.  The  platinum  wire  is  electrically  connected  to  the  HP  3852A 
Data  Acquisition  System.  The  Data  Acquisition  System  will  use  an  internal  ohmmeter  to 
measure  the  resistance  of  the  wire  during  the  procedure.  First  the  line  resistance  of  the 
platinum  wire  is  measured.  The  temperature  of  the  bath  will  be  increased  from  20°C  to 
85°C  with  increments  of  5°C.  When  the  bath  temperature  reaches  steady  state  at  each 
temperature  increment,  the  resistance  of  the  platinum  wire  and  the  bath  temperature  is 
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recorded.  Once  the  data  is  recorded  at  85°C,  the  process  is  repeated  using  temperature 
decrements  of  10°C  until  20°C  is  achieved.  The  bath  temperature  and  platinum  wire 
resistance  are  measure  again  at  each  interval  once  steady  state  is  achieved.  The 
measurements  are  recorded  during  the  temperature  decrements  in  order  to  detect  possible 
hysteresis  in  measurements. 

C.  CALIBRATION  PROCEDURE  IN  THE  MAIN  CHAMBER 

After  the  platinum  wire  is  placed  in  the  main  chamber  and  the  chamber  is  filled 
with  FC-72,  another  calibration  procedure  was  conducted  prior  to  securing  the  heat 
exchanger  to  the  top  of  the  main  chamber.  The  procedure  will  use  the  same  equipment  as 
in  Part  B  with  the  exception  of  the  calibration  bath  and  power  supply.  The  FC-72  will  be 
heated  in  the  main  chamber  by  the  installed  strip  heaters.  The  temperature  of  the  FC-72 
will  start  at  20  °C  and  will  rise  in  5  °C  increments  until  saturation  temperature  of  56  °C  is 
achieved. 

D.  RESULTS  AND  ERROR  ANALYSIS  FOR  CALIBRATION 

The  calibration  data  from  the  test  run  in  the  calibration  bath  is  displayed  in  Table 
4  and  the  plot  shown  in  Figure     .  The  curve  fitting  of  the  acquired  data  resulted  in  the 
following  equation; 

TSUT{=  49.689  •  Rpt wire  -253.836 
Where:  Tsurf     =    Surface  temperature  of  the  platinum  wire  in  °C. 

Rptwire  =    Resistance  of  the  platinum  wire  in  Ohms. 
The  platinum  wire  was  again  calibrated  upon  the  completion  of  the  experiments  to  verity 
proper  operation.  The  resulting  calibration  curve  produced  higher  temperature  values  for 
a  given  wire  resistance.  This  change  in  the  calibration  results  could  be  due  aging  in  the 
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platinum  wire  created  by  the  cyclic  thermal  stresses  imposed  on  the  wire  during  the 
experiments.  The  resulting  calibration  curve  data  is  given  in  Table  7  and  a  graph  of  both 
calibration  curves  is  Figure  18. 
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Resistance 
(Ohms) 

Temperature 
(°C) 

Temperature  1 
(°C) 

Percent 
Error  (%) 

5.533 

21.12 

21.09 

0.1421 

5.6185 

25.35 

25.34 

0.3290 

5.749 

30.46 

31.83 

4.480 

5.812 

34.97 

35.11 

0.3870 

5.924 

40.58 

40.35 

0.5667 

6.027 

45.69 

45.63 

0.1103 

6.129 

50.74 

50.71 

0.0630 

6.240 

56.73 

56.22 

0.8930 

6.342 

61.34 

61.29 

0.0788 

6.437 

66.01 

66.59 

0.8780 

6.536 

70.45 

70.93 

0.6830 

6.636 

75.93 

75.50 

0.5663 

6.721 

80.14 

80.12 

0.0202 

6.835 

85.51 

85.79 

0.3270 

6.643 

76.31 

76.21 

0.1310 

6.427 

65.54 

65.51 

0.0378 

6.234 

55.98 

55.82 

0.0978 

6.006 

44.63 

44.49 

0.3140 

5.818 

35.38 

35.25 

0.3540 

5.615 

25.17 

25.16 

0.0397 

Table  4. 


Calibration  Data  and  Linear  Aggression  Results. 


Resistance 
(Ohms) 

Temperature 
(°C) 

5.631 

21.71 

5.693 

24.90 

5.776 

28.30 

5.834 

31.65 

5.917 

35.931 

5.943 

37.24 

5.968 

38.49 

5.985 

39.145 

6.008 

40.587 

6.072 

43.718 

6.129 

46.277 

6.157 

47.665 

6.168 

48.687 

6.230 

51.77 

6.280 

54.268 

6.289 

54.755 

6.293 

55.019 

Table  5. 


In-Place  Calibration. 


52 


Platinum  Wire  (D=0.05mm)  Temperature  Calibration  Data 
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Figure  15.     Calibration  Curve  of  Platinum  Wire. 
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100 


Calibration  Curves  Fit  Compared  to  Experimental  Data 
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Figure  16.    Calibration  Curve  Fit  Compared  to  Experimental  Data. 
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In-Place  Calibration  Data  Compared  to  Curve  Fit 
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Figure  17.    In-Place  Calibration  Data  Compared  to  Calibration  Curve  Fit. 
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Resistance 
(Ohms) 

Temperature 
(°C) 

5.439 

22.35 

5.590 

30.13 

5.721 

36.48 

5.825 

41.69 

6.00 

50.48 

6.097 

55.33 

6.228 

61.91 

6.289 

65.01 

6.387 

69.98 

6.477 

74.48 

6.565 

79.63 

6.694 

85.11 

6.593 

80.13 

6.399 

70.46 

6.2085 

60.75 

6.01 

51.09 

5.841 

42.13 

5.625 

31.56 

5.423 

21.55 

Table  6. 


Final  Calibration  Data 
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Final  Calibration  Curve 
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Figure  18.        Final  Calibration  Curve. 
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Comparison  of  Calibration  Curves 
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Figure  19. 


Comparison  of  Calibration  Curves. 
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APPENDIX  B.  SAMPLE  CALCULATIONS 

A.         DETERMINATION  OF  EXPERIMENTAL  VALUES 

The  following  sample  calculations  are  the  results  of  the  6th  experiment.  This 
experiment  is  the  boiling  curve  of  FC-72. 

1.  Condenser  temperature: 

The  condenser  temperature  is  the  average  of  four  thermocouples  in  the  region  of 
the  condenser/heat  exchanger. 

(7^  +  7^+7^+70 

1  cond 

4 

Tcond  =  27. 123°C 

2.  Fluid  bulk  temperature: 

The  fluid  temperature  is  the  average  of  four  thermocouple  readings.  The 
thermocouples  are  located  inside  the  main  chamber  approximately  one  inch  below  the 
platinum  wire. 

.   (r4  +  r5+7;+:r7) 

1  bulk  -    ; 

4 

Tbuik  =  53.343  °C 

3.  Current  of  platinum  wire  and  precision  resistor: 

The  current  of  the  platinum  wire  and  the  precision  resistor  is  obtained  from  the 
division  of  the  voltage  drop  across  the  precision  resistor  divided  by  its  resistance  value. 


Ipt-Wire  -  I2Q  - 


_    r  2Q 


R2Cl 


IpLWire  -  I2Q 0.2282 A 

2.03 
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4.  Platinum  wire  resistance  value: 

The  platinum  wire  resistance  is  obtained  from  the  division  of  the  voltage  drop 
across  the  platinum  wire  by  the  platinum  wire  current. 

V 

R_    r  PtWire 
P    ~" 


^PtWire 


■* PtWire 


1.4429 

Kpt.Wire  -  TTTT— 
0.2282 

RptWire  =  6.3234  H 

5.  Platinum  wire  surface  temperature: 

To  find  the  platinum  wire  surface  temperature,  the  calibration  formula  in  the 
APPENDIX  A  is  used  with  the  platinum  wire  resistance  as  an  input. 

TSurf  =  49.689  •  Rpt.wire  -  253.836 
Tsurf  =  (49.689  •  6.3234)  -  253.836 
Tsurf=  60.367  °C 

6.  Heat  flux  from  platinum  wire: 

The  heat  flux  from  the  platinum  wire  is  equal  to  the  division  of  electrical  power 
supplied  to  the  platinum  wire  by  the  platinum  wire  surface  area. 

/-\"  _    2  PtWire       r  PtWire 


A 

"■PtWire 


0.2282*  1.4429 
1.76025  10-5 

2 


Q"  =  18,709.12  W/m 

Where  Apt  wire  =  7t  •  Lptwire  •  DptWire 


Aptwire  =  7C- 0.1018736  -5. 112-  10_6=  1.635-  10-6  m2 
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7.  Piston  displacement: 

The  piston  displacement  is  equal  to  the  cross  sectional  area  multiplied  by  the 
piston  stroke  length.  The  following  calculations  in  parts  7  and  8  are  done  for  the  stroke 
length  of  13.3  mm. 


Piston  —  m  U piston  "  ^Stroke 


Apiston  =   —  •  D*  •  L 


APlston=     ^•(9.525)2«13.3 

Apiston=  947.701  mm3 

8.  Amplitude  of  oscillation  in  the  main  chamber: 

The  amplitude  of  the  oscillation  in  the  main  chamber  is  obtained  from  the  division 
of  piston  displacement  by  the  main  chamber  cross  sectional  area. 

A_    ^Piston 
Amplitude  ~ 

a*b 

Where  a  and  b  are  the  inner  dimensions  of  the  main  chamber 

947.701 
A^puwe      1522984#59377 

Amplitude =  0. 105  mm 

9.  Frequency  calculation: 

The  frequency  of  the  oscillation  is  the  arithmetical  inverse  of  the  oscillation 
period  reading  in  seconds.  The  following  calculation  is  completed  for  a  frequency  of 
1.2  Hz. 

F= X—r^  =  — ! —   =  1.2  Hz 
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B.         DETERMINATION  OF  NUMERICAL  VALUES 

In  this  part,  natural  convection  heat  transfer  around  the  platinum  wire  immersed 
in  FC-72  is  calculated  numerically  by  using  recommended  correlation  by  Kuehn  and 
Goldstein  [Ref.  3]  and  Egger's  results  [Ref.  9]  are  used  for  fluid  properties. 

1.  Film  temperature: 

_    Vbulk   +  ■* PtWire) 


Tfilr 


Tfiln 


2 
(53.343  +  60.367) 


Tfiim  —  56.855  °C 

2.  Thermal  conductivity: 

k_  (0.6033- 0.001 15 «rj 
10.0 

K  =  0.05379  W/m°C 

3.  Liquid  density: 

p,  =  (l  .740  -  0.00261  •  T^  )•  1000 

p,=  1.5916-  103kg/m3 

4.  Kinematic  viscosity: 


v=  1. 203952 -lO^exp 


1058.4109 

.  T,lm  +273.15  , 
v  film  J 


v  =  0.29832-  lO^nvVsec 


5.  Specific  heat: 

Cp  =  (0.241 1 1 1  +  3.70337 •  10"4  •  Tfilm)» 4186 

Cp=  1.09742-  103J/kg°C 
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6.  Thermal  expansion  coefficient: 

3 


0.00261 
(1.740-0.00261.  n  J 


P  =  0.0016398   1/°C 


7.  Thermal  diffusivity: 


a 


P/*C, 


a  =  0.31041  •  10-7m2/sec 


8.  Prandtl  number: 


a 
Pr=  9.6105 


9.  Rayleigh  number: 

_  gr,B»(J?t.Wre-rtott) 


RaD  = 


D 


va 


Pi  Wire 


RaD=  1.502 


10.  Nusselt  number: 


NuD  = 


log 


l  +  i= 


0.518  •Ra^** 


1  + 


"0.559V 


3/5 -.-5/12  V5 


+ 


(o.i./te-r 


NuD=  1.579 
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11.  Heat  flux: 


_  k  •  i\uD  •  [I PlWire     lbulk ) 


D 


Q"=   1 1,640  W/m2 
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APPENDIX  C. 


UNCERTAINTY  ANALYSIS 


A.         UNCERTAINTY  IN  SURFACE  AREA 


Dptwire  =  51.1|im  =  0.0511mm  Wire  Diameter 


Lpt.wtre  =  4.05  lin  =  102.8mm 

coD  =  0.002mm 

col  =  0.06in  =  1.52mm 

Apt.  Wire  =  ft  •  LptWire  •  Dpt.wire 

Apt.wire=^»0. 1028m  •  5.11  •  10"6m 


Aptwire=1.65-  10"5m2 


Wire  Length 

Uncertainty  in  Wire  Diameter 
Uncertainty  in  Wire  Length 
Wire  Surface  Area 


dA 


Pt.Wire 


3D 


=  (x*LPtWlre) 


Pt.Wire 


dA 


Pt.Wire 


dL 


=  (x*E>PtMre) 


Pt.Wire 


Then  the  uncertainty  in  surface  area  is: 

©  A,Pt.Wire  =  t](x  •  LptWlre  •  CO  D  f  +  (tT  •  DptWire  •  CO L  f 


CO 


A. Pt.Wire 


1 Pt.Wire 


[  0.002  y  t  A.524^ 

.0.051lJ      1,102.8 


'■Pt.Wire 


CO 


A,Pt.Wire 


*PtWire 


=  0.04 18  or  4. 18% 


C0A,Pt.wire  =  6.91  ■  10"7  m2 
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B.     UNCERTAINTY  IN  POWER 


co v,2Q  =  0.006%  Uncertainty  in  Precision  Resistor  Voltage  Drop  [Ref.  1 1  ] 

cov,pt.wire  =  0.006%      Uncertainty  in  Platinum  Wire  Voltage  Drop  [Ref.  1 1] 


©r,2q  =  0.02O 


Uncertainty  in  Precision  Resistor 


Q  -  Vpt.Wire  •  IptWire 


Q=^, 


2C2 


Wire 


R 


2Q 


dQ 


dV, 


Pt.Wire 


K2Q 


dV. 


20 


Pt.Wire 
.     ^2Q     J 


dO 


dR 


20 


Y  PtWire       y  2Q 


R 


2Q 


Then  the  uncertainty  in  power  is: 


coQ  = 


"  2fl 
V^2Q 


V   (v. 


CO 


V.Pt.Wire 


V    r 


Pt.Wire 
.     ^2Q 


•CO 


K.20 


+ 


V       *v 

y  Pt.Wire        ¥  2Q 


R 


CO 


20 


2C1 


0), 


Q 


0) 


V, PtWire 
V    'Pt.Wire     J 


v  r 


+ 


6) 


K.2Q 


f«     Y 
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c. 


UNCERTAINTY  IN  HEAT  FLUX 
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Then  the  uncertainty  in  heat  flux  is: 
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=  0.0429  or  4.29% 


D. 


UNCERTAINTY  IN  TEMPERATURE 


The  uncertainty  in  the  thermocouple  temperature  measurements  is  the  sum  of  the 
uncertainties  due  to  the  hardware  and  software,  namely  HP  3  852 A  and  the  software  used 
in  this  study.  According  to  Ref.  1 1,  the  temperature  measurement  accuracy  with  the  HP 
3852A  Data  Acquisition  System  is  ±0.4°C  in  the  temperature  range  of  0  to  300°C. 

cotc  =  0.4°C 
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UNCERTAINTY  IN  WIRE  SURFACE  TEMPERATURE 

As  it  is  given  in  Appendix  A,  Section  C  'Results  and  Error  Analysis'  maximum 


uncertainty  is: 
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^  =  3.031% 


surf 


F. 


©Tsurf  =  0.8696°C        For  the  surface  temperature  of  28.69°C. 
UNCERTAINTY  IN  OSCILLATION  AMPLITUDE 
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G.         UNCERTAINTY  IN  OSCILLATION  FREQUENCY 

Because  of  the  friction  factor  uncertainty  in  the  frequency  is  higher  at  the  lower 
frequencies: 
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APPENDIX  D. 


COMPUTER  PROGRAM 


This  program  will  be  used  for  data  reduction.  It  will  take  the  acquired  data  from 
the  experimental  apparatus  and  calculate  the  heat  flux  and  temperature  overshoot.  The 
heat  flux  as  calculated  by  the  natural  convection  heat  transfer  correlation  will  also  be 
plotted  and  compared  to  the  experimental  results. 


Data  Reduction  Program  During  the  Experiment 

This  part  of  the  program  will  take  the  raw  data,  correct  the  temperatures  using  the 
calibration  equations,  and  send  the  data  into  the  main  part  of  the  program. 
C  =  [-0.0001  1.0058-0.0199; 

-0.0001  1.0072-0.0493; 

-0.0001  1.0068-0.0572 

-0.0001  1.0073-0.0509 

-0.0000  1.0047-0.0028 

-0.0001  1.0070-0.0523 

-0.0001  1.0077-0.0583 

-0.0001  1.0089-0.0797 

-0.0001  1.0073-0.0540 

-0.0001  1.0062-0.0471 

-0.0002  1.0158-0.1899. 

-0.0001  1.1186-0.0809; 

0.0000  1.0000  0.0000; 

0.0000  1.0000  0.0000]; 
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forj=  1:35 

ifG=i) 

load  runOl 
t=run01; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


if  0=2) 

load  run02 
t=run02; 

if  0=3) 

load  run03 
t=run03; 

if  0=4) 
load  nin04 
t=run04; 

if  0=5) 
load  run05 
t=run05; 

if  0=6) 
load  run06 
t=run06; 

if  0=7) 
load  run07 
t=run07; 

if  0=8) 
load  run08 
t=run08; 

if  0=9) 
load  run09 
t=run09; 

if  0=10) 
load  run  10 
t=runlO; 

if  0=H) 
load  run  11 

t=runll; 
if  0=12) 
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end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


load  run  12 
t=runl2; 

if  0=13) 

load  run  13 
t=runl3; 

if  0=14) 

load  run  14 
t=runl4; 

if  0=15) 
load  run  15 
t=runl5; 

if  0=16) 
load  run  16 
t=runl6; 

if  0=17) 
load  run  17 
t=runl7; 

if  0=18) 
load  run  18 

t=runl8; 

if  0=19) 
load  run  19 
t=runl9; 

if  0=20) 
load  run20 
t=run20; 

if  0=21) 
load  run21 
t=run21; 

if  0=22) 
load  run22 
t=run22; 

if  0=23) 
load  run23 
t=run23; 
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end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


end; 


if  0=24) 

load  run24 
t=run24; 

if  0=25) 
load  run25 
t=run25; 

if  0=26) 
load  run26 
t=run26; 

if  0=27) 
load  run27 
t=run27; 

if  0=28) 
load  run28 
t=run28; 

if  0=29) 
load  run29 
t=run29; 

if  0=30) 
load  run30 
t=run30; 

if  0=31) 
load  run3 1 
t=run31; 

if  0=32) 
load  run32 
t=run32; 

if  0=33) 
load  run33 
t=run33; 

if  0=34) 
load  run34 
t=run34; 
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if  0=35) 
load  run35 
t=run35; 
end; 


fori=  1:14 

T(iJ)  =  C(i,l)*t(ir2  +  C(i,2)*t(i)  +  C(i,3); 
end; 


end; 


This  part  of  the  program  will  be  used  in  data  reduction  for  the  experiment. 

First  the  data  is  loaded  into  the  active  part  of  the  program. 

Q  =  zeros(l,35); 
Tsh  =  Q; 

forj  =  l:35 
The  first  calculation  will  be  to  determine  the  temperature  of  the  condenser. 

Temp  =  T(:  j); 

Tc  =  mean(Temp(l:4)); 
Next  the  temperature  of  the  bulk  fluid  will  be  calculated. 

Tb(j)  =  mean(Temp(5:7)); 

The  voltage  drop  will  be  used  to  find  Resistance  of  Platinum  wire. 
Vr(j)  =  Temp(13); 
Vp(j)  =  Temp(14); 
Ic(j)  =  Vr(j)/2.03; 

ip(j)=ic0); 
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Rp(J)  =  vPG)/ipG); 

The  resistance  of  the  platinum  wire  will  be  used  to  calculate  the  actual  surface 
temperature  of  the  platinum  wire. 

TsurfG)  =  49.989*RpG)  -  253.847; 

The  film  temperature  of  the  fluid  in  direct  contact  with  the  wire  is  calculated  using  the 
surface  temperature  and  the  bulk  temperature. 

TfilmG)  =  (Tsurf(j)  +  TbG))/2; 

The  wall  superheat  needs  to  be  calculated.  This  value  is  a  measure  of  the  temperature 
difference  between  the  platinum  wire  and  the  saturation  temperature  of  the  FC-72. 

TshG)  =  TsurfG)  -56.1; 

The  heat  flux  generated  by  the  wire  is  calculated  based  on  the  power  output  of  the  wire 
and  its  area. 

QG)  =  (IpG)*VpG))/l-76025e-5; 

The  two  values  that  need  to  be  saved  from  each  data  run  are  the  heat  flux  Q  and  the  wall 
superheat  Tsh. 

This  part  of  the  program  will  calculate  the  natural  convection  heat  flux 
based  on  the  Kuehn  and  Goldstein  correlation. 

ifj<=8 

Thermal  conductivity  calculation.  (W/cmA2) 

k  =  ((0.6033  -  0.001 15*Tfilm(j))/10); 

Liquid  density  calculation.    (kg/mA3) 

p  =  (1.74  -  0.00261*TfilmG))*1000; 

Kinematic  viscosity  calculation.    (mA2/sec) 

v  =  ((1.203952e-08)*(exp(1058.4109/(TfilmG)+273)))); 

Specific  heat  calculation.     (J/kg  C) 

Cp  =  (0.241 111  +  3.70337e-04*TfilmG))*4186; 
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Thermal  expansion  coefficient.    (1/  C) 

B  =  0.00261/(1.740- 0.0025  l*Tfilm(j)); 
Thermal  diffiisivity. 

al  =  k/(p*Cp); 
Prandlt  number. 

Pr  =  v/al; 

Rayleigh  number. 

RaG)  =  (9.817*B*(Tsurf(j)  -  Tb(j))*(0.05e-03)A3)/(v*al); 

Nusselt  number. 

NuG)  =  2/log(l+(2/((0.518*RaG)A0.25*(l+(0.559/Pr)A0.6)A-0.4167)A15 

+  (0.1*RaG)A0.333)A15)A0.0667)); 


Heat  flux. 

qG)  =  k*NuG)*(TsurfG)  -  TbG))/0.05e-03; 

else 

zG)  =  0; 

end 


end; 

This  part  of  the  program  will  reorganize  the  data  and  prepare  it  to  be  graphed  and 
displayed  so  the  values  can  be  compared  to  one  another. 


KG  =  Tsh(l:8); 
KGQ  =  Q(1:8); 


Tbl  =  Tb(l:22);  Tb2  =  Tb(22:35); 
Tl  =  Tsh(l  :22);  T2  =  Tsh(22:35); 
Q1=Q(1:22);        Q2  =  Q(22:35); 
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format  short  e 

semilogy(T  1,Q1,V,T  1  ,Q  1  /-•,T2,Q2,'o,,T2,Q2,'- 

,,KG,q,V,KG,q,,:,,Tbl,Ql,,+,,Tb2,Q2,,,),grid; 

axis([0,90, 1 0000, 1 000000]); 

title('Heat  flux  vs.  Wall  superheat'); 

xlabel('Wall  Superheat  deg  C); 

ylabel(*Wire  Heat  Flux  W/cmA2'); 

This  part  of  the  program  will  display  the  two  heat  flux  values  and  the 
temperature  overshoot  value. 

Final  =  ([Q^Tsh1]) 
Flux  =  ([KGQ'^Cq']) 

This  will  save  the  graph  to  a  file  for  easy  printout  later. 

print  -dps2  plotl 
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APPENDIX  E.  EXPERIMENTAL  DATA 

A.         EXPERIMENTAL  DATA 

1.  The  results  of  each  experimental  run  are  displayed  in  Table  6.  Most 
experimental  runs,  which  were  not  completed,  were  halted  due  to  leaks  developing  in  the 
main  chamber.  The  highly  wetting  FC-72  was  able  to  seep  through  the  seams  of  the  main 
and  front  chambers.  Re-applying  silicon  sealant  after  each  set  of  experiments  to  seams  of 
main  chamber  halted  the  leaks. 

2.  Figures  20-32  are  the  heat  flux  vs.  wall  superheat  plots  from  the 
experimental  runs  that  were  not  displayed  in  Chapter  IV. 
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Experiment 

Results 

1 

Stopped  --leaks 

2 

Stopped  -  leaks 

3 

Stopped  -  broken  wire 
leading  to  power  supply 

4 

Completed 

5 

Stopped  -  leaks 

6 

Completed 

7 

Stopped  --pump failed 

8 

Stopped  -  leaks 

9 

Stopped  -  leak  in  heat 
exchanger 

10 

Stopped  -  leaks 

11 

Stopped  -  incorrect 
piston  stroke  length 

12 

Completed 

13 

Completed 

14 

Completed 

15 

Completed 

16 

Completed 

17 

Completed 

18 

Completed 

19 

Completed 

20 

Completed 

21 

Completed 

22 

Completed 

23 

Completed 

24 

Completed 

25 

Completed 

26 

Completed 

27 

Completed 

28 

Completed 

29 

Completed 

30 

Stopped  -  leaks 

31 

Stopped  -  leaks 

32 

Stopped  -  leaks 

Table  7. 


Status  of  Experimental  Runs. 
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Heat  flux  vs.  Wall  superheat 
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Figure  20.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 1 1 7mm  and  a  Frequency  of  1 .0Hz. 
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Figure  2 1 .        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.11 7mm  and  a  Frequency  of  1 .2Hz. 
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Heat  flux  vs.  Wall  superheat 
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Figure  22.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.105mm  and  a  Frequency  of  0.6Hz. 
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Figure  23 .        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.105mm  and  a  Frequency  of  1.0Hz. 
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Figure  24.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 1 05mm  and  a  Frequency  of  1 .2Hz. 
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Figure  25 .        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.125mm  and  a  Frequency  of  0.6Hz. 
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Figure  26.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.125mm  and  a  Frequency  of  1 .0Hz. 
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Figure  27.         Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 125mm  and  a  Frequency  of  1 .2Hz. 
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Figure  28.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 1 1 7mm  and  a  Frequency  of  1 .0Hz. 
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Figure  29.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 1 1 7mm  and  a  Frequency  of  1 .2Hz. 
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Figure  30.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0. 1 1 7mm  and  a  Frequency  of  0. 1 1 7mm. 
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Figure  3 1 .        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.252mm  and  a  Frequency  of  0.3Hz. 
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Figure  32.        Boiling  Curve  of  FC-72  with  an  Oscillation  of  Amplitude 
0.252mm  and  a  Frequency  of  1.5Hz. 
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